Current self-oscillations, spikes and crossover between charge monopole and dipole 

waves in semiconductor superlattices 
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Self-sustained current oscillations in weakly-coupled superlattices are studied by means of a 
self-consistent microscopic model of sequential tunneling including boundary conditions naturally. 
Well-to-well hopping and recycling of charge monopole domain walls produce current spikes -high 
frequency modulation- superimposed on the oscillation. For highly doped injecting contacts, the self- 
oscillations are due to dynamics of monopoles. As the contact doping decreases, a lower- frequency 
oscillatory mode due to recycling and motion of charge dipoles is predicted. For low contact doping, 
this mode dominates and monopole oscillations disappear. At intermediate doping, both oscillation 
modes coexist as stable solutions and hysteresis between them is possible. 
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Solid state electronic devices presenting negative dif- 
ferential conductance, such as resonant tunneling diodes, 
Gunn diodes or Josephson junctions are nonlinear 
dynamical systems with many degrees of freedom. They 
display typical nonlinear phenomena such as multista- 
bility, oscillations, pattern formation or bifurcation to 
chaos. In particular, vertical transport in weakly cou- 
pled semiconductor doped superlattices (SLs) has been 
shown to exhibit electric field domain formation P-I|, 
multistability ||, self-sustained current oscillations |3-3|, 
and driven and undriven chaos [js|— . Stationary electric 
field domains appear in voltage biased SLs if the doping 
is large enough ||. When the carrier density is below 
a critical value, self-sustained oscillations of the current 
may appear. They are due to the dynamics of the domain 
wall (which is a charge monopole accumulation layer or, 
briefly, a monopole) separating the electric field domains. 
This domain wall moves through the structure and is pe- 
riodically recycled. The frequencies of the corresponding 
oscillation depend on the applied bias and range from 
the kHz to the GHz regime. Self-oscillations persist even 
at room temperature, which makes these devices promis- 
ing candidates for microwave generation Q. Theoreti- 
cal and experimental work on these systems have gone 
hand in hand. Thus the paramount role of monopole 
dynamics has been demonstrated by theory and exper- 
iments. Monopole motion and recycling can be experi- 
mentally shown by counting the spikes -high frequency 
modulation- superimposed on one period of the current 
self-oscillations: current spikes correspond to well-to- well 
hopping of a domain wall through the SL. In typical ex- 
periments the number of spikes per oscillation period 
is clearly less than the number of SL wells |Q,^2|. It 
is known that monopoles are nucleated well inside the 
SL |?],|| so that the number of spikes tells over which 
part of the SL they move. Other possible waves, such 
as the charge dipole waves appearing in the well-known 



Gunn effect, are nucleated at the emitter contact |L3| . 
Had they been mediating the self-oscillation, the number 
of current spikes would be comparable that of SL wells. 

In this letter we study the non-linear dynamics of SLs 
by numerically simulating the model proposed in Ref. 
|14| . Our simulations show self-sustained oscillations of 
the current and current spikes reflecting the motion of 
the domain wall as observed experimentally. Further- 
more, when contact doping is diminished, we predict a 
crossover from monopole to dipole self-oscillations resem- 
bling those in the Gunn effect |13|. Indeed, our results 
show for first time that there is an intermediate range of 
contact doping and a certain interval of external dc volt- 
age for which monopole and dipole self-oscillations with 
different frequencies are both stable. Hysteretic phenom- 
ena then exist. 

1. Model and superlattice sample. Our self-consistent 
microscopic model of sequential tunneling includes a de- 
tailed electrostatic description of the contact regions and 
SL . It consists of a system of 3 A + 8 equations for 
the Fermi energies and potential drops at the N wells, 
the potential drops at the barriers and at the emitter 
and contact layers, width thereof, charge at the emit- 
ter and total current density. These equations comprise 
Ampere current density balance and Poisson equations, 
conservation of the global charge, and the overall volt- 
age bias condition. Dynamics enters the model through 
Ampere's law for the total current density J = J(t), 



J — Ji-l. 



d dt 



(1) 



which is equivalent to a local charge continuity equa- 
tion [pr[ . Here Ji— i,t is the tunneling current density 
through the ith barrier of thickness d, evaluated by using 
the Transfer Hamiltonian approach The last term 

in (|l]) is the displacement current at the ith barrier where 
the potential drop is Vi and e is the static permittivity. 
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Our numerical simulations (of the 37V+8 coupled equa- 
tions) have been performed for a 13.3 nm GaAs/2.7 nm 
AlAs SL at zero temperature consisting of 50 wells and 
51 barriers, as described in [12). Doping in the wells 
and in the contacts are N w = 2 x 10 10 cm -2 and 
N c = 2 x 10 16 cm~ 3 respectively. Notice that the typical 
experimental value is N c = 2 x 10 18 cm~ 3 0,|l2|. For this 
value, we find current self-oscillations due to monopole 
dynamics with very small superimposed current spikes. 
Since the origin of such spikes is the same as for smaller 
N c (for which spikes are larger and bistability of oscilla- 
tions is possible), we choose not to present data compa- 
rable to experiments in this paper (see Ref. |l2| for the 
relevant experimental data). 

2. Monopole- mediated self-oscillations of the current. 
Fig. |](a) depicts the current as a function of time for 
a dc bias voltage of 5.5 V on the second plateau of the 
SL I — V characteristic curve. J(t) oscillates periodically 
at 20 MHz. Between each two peaks of J(t), we observe 
18 additional spikes. The electric field profile is plotted 
in Fig. ||(b) at the four different times of one oscillation 
period marked in Fig. |l|(a). There are two domains of 
almost constant electric field separated by a moving do- 
main wall of (monopole) charge accumulation (which is 
extended over a few wells) . Monopole recycling and mo- 
tion occur on a limited region of the SL (between the 
30th and the 50th well) and accompany the current os- 
cillation . Well-to- well hopping of the domain wall is 
reflected by the current spikes until it reaches the 46th 
well which is close to the collector. Then the strong in- 
fluence of the contact causes that no additional spikes 
appear. Instead the current rises sharply triggering the 
formation of a new monopole closer to the emitter con- 
tact but well inside the SL; see Figs. [j](a) and (b). The 
number of wells traversed by the domain wall (almost) 
coincides with the number of spikes per oscillationpc- 
riod, a feature not found in previous models. Fig. |l)(b) 
shows the recycling of a monopole: between times (1) 
and (3) there is a single monopole propagating towards 
the collector; at (4) a new monopole is generated at the 
middle of the structure and the old one collapses at the 
collector. It is interesting to realize that the region near 
the emitter does not have a constant electric field profile 
due to the large doping there (its Fermi level is well above 
the first resonant level of the first well). This produces a 
large accumulation layer. 

3. Current spikes. What is remarkable in Fig. |l](a) (as 
compared to previous studies) are the spikes superim- 
posed near the minima of the current oscillations. Such 
spikes have been observed experimentally and attributed 
to well-to- well hopping of the domain wall Jl2|,[l5| . They 
are a cornerstone to interpret the experimental results 
and in fact support the theoretical picture of monopole 
recycling in part (about 40%) of the SL during self- 
oscillations. The identification between number of spikes 
and of wells traversed by the monopole rests on volt- 
age turn-on measurements supported by numerical sim- 
ulations of simple models during early stages of station- 



ary domain formation |15|]. These models do not pre- 
dict spikes superimposed on current self-oscillations due 
to monopole motion j|j7]|l(|. To predict large spikes, a 
time delay in the tunneling current |1^] or random dop- 
ing in the wells [|l7]] have to be added. Unlike these mod- 
els, ours reproduces and explains spikes naturally thereby 
supporting their use to interpret experimental results. 

Fig. ||(a) depicts a zoom of the spikes in Fig. |l|(a) . They 
have a frequency of about 500 MHz and an amplitude of 
2.5 /j,A. Fig. ||(b) shows the charge density profile at four 
different times of a current spike marked in Fig. ||(a). No- 
tice that the electron density in Fig. ||(b) is larger than 
the well doping at only three wells (40, 41 and 42) during 
the times recorded in Fig. |^(a) . The maximum of electron 
density moves from well 40 to well 41 during this time 
interval so that: (i) tunneling through the 41st barrier 
(between wells 40 and 41) dominates when the total cur- 
rent density is increasing, whereas (ii) tunneling through 
barriers 41 and 42 is important when J(t) decreases. The 
contributions of tunneling and displacement currents to 
J(t) in Eq. (|l]) are depicted in Figures ||(c) and (d). 

More generally, the spikes reflect the two-stage hop- 
ping motion -fast time scale- of the domain wall: at 
time (1) (minimum of the current), the charge accumu- 
lates mainly at the i-th well. As time elapses, electrons 
tunnel from this well to the next one, the (i+l)-st, where 
most of the charge is located at time (3) (maximum of the 
current). This corresponds to a hop of the monopole. As 
the monopole moves, it leaves a lower potential drop on 
its wake. The reason is that the electrostatic field at the 
(i+l)-st well and barrier become abruptly flat between 
times (1) and (3), as they pass from the high to the low 
field domain. This means that a negative displacement 
current has its peak at the (i+1 )-st barrier, near the wells 
where most of the charge is. Between times (1) and (3), 
the tunneling current is maximal where the displacement 
current is minimal and the total current increases. After 
that, some charge flows to the next well [time (4)] but 
both, tunneling and displacement currents, are smaller 
than previously. This occurs because the potential drop 
at barrier (i+2) (in the high field domain) is larger than 
at barrier (i+1). Then there is a smaller overlap be- 
tween the resonant levels of nearby wells -the tunneling 
current decreases - and the displacement current and, 
eventually, J(t) decreases. This stage lasts until well i is 
drained, and most of the charge is concentrated at wells 
(i+1) (the local maximum of charge) and (i+2) (slightly 
smaller charge). Then the next current spike starts. 
4- Dipole self-oscillations of the current. An advantage 
of our present model over other discrete ones J|,[l6|,[l8) is 
our microscopic modeling of boundary conditions at the 
contact regions. Thus we can study what happens when 
contact doping is changed. The result is that there ap- 
pear dipole-mediated self-oscillations as the emitter dop- 
ing is lowered below a certain value. There is a range of 
voltages for which dipole and monopole oscillations co- 
exist as stable solutions. This range changes for different 
plateaus. When the emitter doping is further lowered, 
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only the dipole self-oscillations remain. FigJ| presents 
data in the crossover range (below N c = 4.1 x 10 16 cm -3 
and above N c = 1.7 x 10 16 cm~ 3 for the second plateau), 
for the same sample, doping and bias as in Figs. ^ and |^. 
Except for the presence of spikes of the current, dipole re- 
cycling and motion in SLs are similar to those observed in 
models of the Gunn effect in bulk GaAs |b| . These self- 
oscillations have not been observed so far in experiments 
due to the high values of the contact doping adopted in 
all the present experimental settings. Notice that cur- 
rent spikes appear differently than in the monopole case, 
Fig. |l|(a). The main difference is that now there are 
many more current spikes, 36, for the dipoles recycle at 
the emitter and traverse the whole SL. See Figs. ||(b) 
and (c). Charge transfer and balance between tunneling 
and displacement current during a spike are similar to 
those occurring in monopole oscillations. For a simpler 
model the velocity of a charge accumulation layer 
(belonging to a monopole or a dipole) has been shown to 
approximately obey an equal area rule. Then monopole 
and dipole velocities are similar but a monopole traverses 
a smaller part of the SL than a dipole does. Therefore 
dipole oscillations have a lower frequency than monopole 
ones. Our results agree with this: the frequency of the 
dipole oscillations discussed above is about 8 MHz, 40% 
the frequency of monopole oscillations. 

Dipole self-oscillations have also been predicted to oc- 
cur in weakly-coupled SLs as the result of assuming a lin- 
ear current - field relation at the injecting contact on a 
simpler model [[l6 19 1. Since such ad hoc boundary condi- 
tion has no clear relation to contact doping, no crossover 
between different oscillation types could appear in that 
work. 

5. Multistability. Monopole and dipole waves coexist 
in both the first and the second plateaus. The time- 
averaged current as a function of dc voltage in the first 
plateau (whose crossover range is below N c — 2.1 x 
10 16 cm~ 3 and above N c = 1.5 x 10 16 cm~ 3 ) has been plot- 
ted in Fig. H Notice that the average current of dipole 
oscillations is lower than that of monopole oscillations. 
Previous studies for Gunn oscillations (p| found that 
large dipole waves appear only for small current values, 
whereas monopole recycling requires current values near 
the maximum of the current-field characteristic curve. 
Let us start at a bias of 0.5 V (for which the stationary 
state is stable) and adiabatically increase the voltage. 
The result is that we go smoothly from the stationary 
state to the fast monopole self-oscillation at about 1.3 V. 
This branch of oscillatory states eventually disappears at 
about 2.6 V. If we now adiabatically lower the bias, we 
reach a slow dipole self-oscillation at about 2.4 V. There 
is a small hysteresis loop between dipole oscillations and 
the stationary state between 2.4 V and 2.6 V: the former 
may start as a subcritical Hopf bifurcation. About 0.8 V 
the dipole oscillation disappears and we are back at the 
stable stationary state. We therefore find the hysteresis 
loops marked by arrows in Fig. ^. 

In conclusion, we have dealt with self-sustained os- 



cillations of the current in SLs whose main mechanism 
is sequential tunneling. Depending on contact doping, 
these oscillations may be due to recycling and motion 
of two different charge density waves: monopoles and 
dipoles. Experimentally, only the monopole oscillations 
have been observed, for the contacts doping is usually 
set to values which are too high. The dipole-like oscilla- 
tions could be observed constructing samples with lower 
doping at the contacts. In fact, as the doping of the con- 
tacts is reduced, we predict current oscillations due to 
dipole charge waves. The crossover between both types 
of self-oscillations occurs at intermediate emitter doping 
values for which stable monopole and dipole oscillations 
coexist. Then the diagram of average current versus 
voltage is multivaluated, presenting hysteresis cycles and 
multistability between monopole and dipole oscillations 
(and between oscillatory and stationary states). The 
time-resolved current in the oscillatory modes presents 
a number of sharp spikes. They occur because well- 
to-well hopping of charge accumulation layers occurs in 
two stages: during the stage where the current rises, 
charge is mainly transferred through a single barrier. 
The charge is transferred through two adjacent barriers 
at the stage in which the current decreases. All these 
properties form the basis for possible applications of SLs 
working as multifrequency oscillators in a wide range of 
frequencies. Quantitative description of such multifre- 
quency oscillators requires calculation of typical output 
power characteristics and noise levels. This is the pur- 
pose of a future work. 
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FIG. 1. (a) Self-sustained oscillations of the total current 
through the SL due to monopole recycling and motion. Bias 
is 5.5V and emitter doping, N c = 2 x 10 16 cm -3 , (b) Electric 
field profiles at the times marked in (a) during one period of 
the current oscillation. 
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FIG. 2. (a) Zoom of Fig. [l] showing the spikes of the cur- 
rent, (b) Electron density profiles (in units of the doping at 
the wells), (c) tunneling current, and (d) displacement current 
within the monopole at the times marked in (a). 
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FIG. 3. (a) Dipole-mediated self-oscillations of the current 
at 5.5V" for iV c = 2 x 10 16 cm" 3 . (b) Detail of the current 
spikes, (c) Electric field profiles at the times marked in (a). 
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FIG. 4. I-V characteristics at the first plateau, for both 
sweep directions showing bistability between self-oscillations 
mediated by monopole and by dipole dynamics. Notice the 
hysteresis cycle. 
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